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ABFTRACT 


Analytical and cxpentacntal studies have 
been undertaken on propellers operating in 
the unsteady flow field produced by inter- 
action effects due to the fuselage, wing, and 
nacelles. Methods have been developed and 
verified experitacntally for determining the 
velocity field in which a propeller operates 
as well as its aerodynamic and dynamic re- 
sponse to this unstcadv environment. Methods 
are presented for predicting the net thrust 
of a propeller-wing-bodv combination as well 
as the unsteady tnrust and torque acting on 
the propeller. Sample calculations as well 
as wind lunnel and flight test results are 
presented which illustrates tne sensitivity 
of a propeller to the flow field in which it 
is operating. 



THE AERODYNAJJIC AND DYf.’AMlC BEIUVIOR of pro- 
pellers for general aviation aircraft are of 
inportance to the operation of tnese aircraft. 
A propeller vhlch is poorly matched to its 
flow field or engine can be inefficient anc 
subject to excessive vibratory stresses. The 
study of propeller aerodynamics and dynamics 
described here is divided into three parts. 

The first part deals with predicting the flow 
field in which the propeller is to operate. 
Here, a numerical method is described for 
calculating the velocity vector at any radial 
and azimuthal location in the propeller plane 
as a function of wing-fuseiage-nacclle geo- 
metry. The second part treats botn the 
steady and unsteady airloads produced by the 
propeller blades moving through the spatially 
varying velocity field. The third, and final 
section, examines the structural dynamics of 
the blades and presents methods for predict- 
ing their normal modes. In essence this 
paper is a brief summary of references 1, 2, 
and 3. For more details, a study of these 
references is recommended. 

PROGRAM FOR PREDICTION OF PROPELLER FLOW FIELD 


Potential flow methods, at least for 
tractor propellers, will accurately predict 
the velocity field in which a propeller 
operates since the effect of viscosity ahead 
of a body is usually negligible. However, 
even with the simplifications afforded by 
potential flow, the calculations for a 
fuselagc-wing conbination can prove tedious. 
Therefore, part of this study investigated 
the accuracy with which one needs to model 
the complete aircraft geometry in order to 
obtain a sufficiently precise velocity field 
at the propeller plane. 
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M"'.ERICAL MODEL 


^ Following the lead of referencesv^ 
( 5 } the fuselage surface is divided into 


0^"an^ 






/ number of small panels as illustrated in a 
K general sense in figure 1. Each panel is 

covered with a constant distributed source 
strength per unit area denoted by o(I). The 
velocity induced by the panel on the 
surface of that panil directed outwa’-d and 
normal to the panel at tlie control point will 
be given by 



V .2SI1 

K 2 


( 1 ) 


th 

For anv other panel, sav the J panel, the 
total source strength over the panel is given 
by 

Q(J) - o(J) S(J) (2) 

In order to calculate the velocity induced at 
the panel by Q(J), Q(J) is taKen to be a 
point source located on the panel control 
point. Letting (})(J) denote the velocity 
potential associated with Q(J) and K(I), the 
unit vector nornal to the It-h panel and 
directed outward, the velocity induced by 
Q(J) nornally outward at the panel will 
be 

Vjj(I.J) - grad ^(J) » N(I) (3) 

If V is the free-stream velocity vector and 
V^,(I) is the wing Induced velocity vector at 
the panel, the conponents of these 
vectors normally outward at the panel 
control point will be 

V • N(I) 

and, 

V^,CI) • N(I) 

w 


The nornal velocity must vanish at panel 1 
if it is a solid boundary or must equal the 
specified through velocity, q, normal across 
tne panel if it is a relaxed boundary which 
models a region of through flow. Satisfying 
the appropriate boundary condition on panel 
1, it follows that 


V • N(I) + V (I) • N(I) + + i V (I,J) 

J=1 


0 (solid boundary) .j , j. 
q (relaxed boundary) 


For a point source of strength Q, .‘he 
velocity potential ia given by. 




Q 


(M 


( 5 ) 


B. 


where r is the radial distance from the 
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Writing, 


tj 


N(I) 


V iu + jv + kw 

(a) 

lu d) + 3 V,.(I) + iw..d) 

V V V 

(b) (6) 

in (I) 4 kn (I) + to,d) 

X y z 

(c) 


Equation 4 can be written as, 

|u + u (I)]n (I) + (v + v^.(I))n^.(l) 

" '' ' (7) 

+ [w + w^(I)]n^(I) + 

1 I £(J)_sm (ixd) - x(J)]n (I) 

+ ly(i) - y(J))ny(i) + lz(i) 

/ M /T\1 J 0 (solid boundary) 

- z(J)]n^(I)i “ * q (relaxed boundary) 


where, 

R(I,J) = {[xd) - x(J)]^ + Iy(i) - y(J)] 

1 

2 2 

+ [z(D - z(J)] } 

Letting I =• 1. 2, ..., K, the foregoing 
represents a set of ^ simultaneous equations 
for the unknowns o(I). In order to solve 
this set of equations, the surface of a given 
body is divided into N panels that can be 
either four-sided or triangular. The •’ask of 
paneling is a formidable one requiring care 
in the indexing of the corners. In the case 
of the non-planar, four-sided panel, the unit 
normal is obtained from the vector product of 
the diagonals. The panel area is then taken 
as the area of the panel projected on a plane 
perpendicular to this unit vector. Arithmetic B 
averages of the x, y, and z positions of the 
four corners define the location of the con- 
trol point and concentrated source for this 
type of panel. Control points and concen- 
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trated sources are located at the centroid of 
triangular panels. 

Once the 0(1) values are obtained, the 
velocity field at any point can be deternined 
by combining the free stream velocity and 
local wing induced velocity vectors with the 
gradient of 4> obtained from equations 2 and 5. 

The total 4), of course, is obtained by sum- 
ming over the index I. The complete program 
which has been developed will nandle both 
wing-fuselage and wing-fuselage-nacelle com- 
binations as shown in figure 2. By employing 
the relaxed boundary condition in equation 7 
on particular panels, it also allows for flow 
through a panel to represent air intakes or 
exhausts. The lifting wing is modeled simply 
by a single horseshoe vortex placed along the 
quarter-chord line (unswspt) and trailing from 
the ti/ 4 spanwise station. This simple model 
is considerably easier to handle computation- 
ally compared to vortex lattice models and 
results in upwash velocities at typical pro- 
peller locations that are within a few percent 
of more exact calculations. As equation 7 
indicates, the effect of the vortex on the 
body is considered in determining the o(I) 
values. However, the effect of the fuselage 
or nacelles on the wing is not taken into 
account. 

The final program has been checked 
against closed-form solutions and the results 
compare favorably, at least for simple, non- 
lifting body shapes. It was then applied to 
typical fuselage shapes with a certain amount 
confidence. These results will be described 
later. 

METHOD FOR PREDICTING AERODWAMIC LOADS ON A 
PROPELLER 

A propeller operating at an angle of 
attack or in the presence of a fuselage 
experiences tine-dependent acrod>namlc leads 
as its blades pass through a spatially vary- 
ing velocity field. Generally, both the 
magnitude and direction of the velocity field 
will be a function of blade position. Thus 
blade section lift coefficients and resultant 
velocities vary periodically with the blade 
azimuth angle. B 

In order to predict the unsteady forces 
on a propeller blade, it is convenient to 
consider first the steady case. A section 
of a propeller blade as viewed looking toward 
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the hub along the blade is pictured in figure 
3. The blade section is operating under the 
influence of three velocity conponents: the 

free stream velocity, V; the linear velocity 
due to rotation, vr; and the induced velocity, 
w. The modem trend is to calculate the in- 
duced velocity w by applying the Biot-Savart 
relationship to the propeller's helical 
trailing vortex system. Even here, however, 
it is usually necessary to assume a shape to 
the wake before numerically integrating the 
Biot-Savart equations. A free-wake analysis, 
where the location of the trailing vortex 
system and its induced velocities are mutually 
dependent, requires a prohibitive amount of 
computer time. 

A propeller application of classical vor- 
tex theory incorporating corrections for pro- 
file drag, finite section thickness, and finite 
chord appears to predict propeller thrust and 
torque in close agreement with experimental 

measurements over a wide range of propeller /' 

geometries. This method, used here, is much 
simpler to apply by comparison to vortex lat- 
tice methods. The details can be found in 
references 1, 6, and 7 and will not be repeated 
here. However, a few points will be noted. 

1. An accurate empirical correlation of 
section C^j with a and Mach number for Clark-Y, 
series 16 and NACA 4-digit airfoil sections is 
to be found in reference 1. 

2. The tangential component of w is 
related to the bound circulation through 
Prandtl’s tip loss factor F. Reference 1 
shows that this closed form approximation to 
Goldstein's kappa factor gives predicted 
results very close to those obtained using 
kappa. 

3. Normality between w and the resultant 
velocity is assumed in the propeller plane. 

This assumption is justified in reference 7. 

4. The tip angle of the helical vortex 
wake necessary to the calculation of F is 
obtained by iterating on w. Figures 4 and 5 
are two typical results of steady propeller 
performance predictions to be found in refer- 
ence 1. The agreement shown in figure 5 
between predictions and experimental results 
IS particularly impressive considering the 

fact that the resultant tip Mach number ex- B. W. McCormick 

ceeds 1.0 for advance ratios greater than 

1 . 2 . 
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CALCULATION OF UNSTEADY PROPELLER AERODYNAMIC 
LOADS 


The prediction of unsteady aerodynamic 
loads builds on the steady case. Generally, 
under the Influence of angle of attack, wing 
upwash,and body interference, the propeller 
operates in a velocity field which varies both 
radially and circumferentially. The approach 
taken here is to express the axial and tan- 
gential velocity components relative to the 
propeller at a given radius as the sum of 
harmonics of i|/, the azimuth angle, is de- 
fined in figure 6. 

• Let and denote the K harmonics 
of the axial and tangential velocity components 
respectively at a particular radius. That is, 

'^AN ' ^^AN ^AH 

( 8 ) 

Vtjj *= cos Ntj; + b^ sin (b) 

Combining and Vjjj vectorially, the N^^ 

harmonic or the velocity vector normal to the 
steady resultant velocity Ve, shown in figure 
7, can be written as, 

''nn “ ''an W + sin + o^) (9) 


Combining 


V /V 

nn''^ 


+ 


+ 


equations 8 and 9 
ISaj. cos «) + a^) 

sin (({. + a^)] 
f^AN “i^ 

sin Hil> 


gives. 


cos Ni/^ 


- sin 


(<f 


ao) 


By comparison to an airfoil undergoing 
a pure heaving motion, Vnk corresponds to the 
heaving velocity h. The heaving acceleration 
h ir given by differentiating V..,, with respect 
to time. Since to = d,^/dt it follows that, 

ii/V = -N'jj[a^,, cos « + 

- a„, sin (C + o^)] sin Ni< (11) 

Ti\ X 

+ Ntj{b,„ cos << + Q.) 

A;> 1 
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- bjj; sin (<> + a^)J cos Kip 



Knowing h and h, one can calculate the un- 
steady lift on the heaving airfoil from 


L = h + TOVcC(k)h (12) 

The above can be found in reference 8. C(k) 

IS a complex function known as Tneordorsen' s 
function while k is the so-calied reduced 
frequency defined b>, 

k “ ^ a3) 

Here o) is the circular frequency of the har- 
monic heaving notion which in this case is 
equal to the rotational speed of the pro- 
peller. The complex function C(k) is fre- 
quently seen expressed as, 

C(k) = F(k) + lG(k) 0.U) 

C(k) is presented in figure 8. 

Thus, 

2 

L/h = irpVcF + i[ ^ — + TipVcG] (15) 

This can be written in the dimensionless form. 


L 

TipVch 

where: 


/ + (j+ G)-^ 


i<P 


4 , = 


-1 

tan 



(16) 


Thus, equation 16 relates the amplitude and 
phase of che unsteady section left, L, to 
the heaving velocity h. <fi is the angle by 
which the lift force leads h. Substituting 
(equation 10) for h thus leads to the 
amplitude and phase for the sectional un- 
steady lift on the propeller blade. This lift 
can be resolved into thrust and torque com- 
ponents and then intcgtcfed over the blade 
length to give the total unsteady forces and 
moments acting on the blade. B. 

The computer program which accomplishes 
the foregoing is described in reference 1. 

As input this program requires the axial and 7 

ta.igentlal velocity components relative to the 
propeller axis and in its disk plane. These 
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are obtained directly as output fron the pro- 
graTO developed in reference 3. These are 
then Fourier-analyzed by the former program 
to obtain the a^, Ajfj, b^jj, and by^ coef- 
ficients. Generally, the first iiarnonic is 
the predominant one. In the results to follow, 
only the first two harmonics are used although 
the program is capable of handling up to eight. 

RESULTS OF UNSTEADY AIRLOAD CALCULATIONS 

A propeller m a uniform flow inclined at 
an angle to the flow of a experiences only a 
first harmonic variation in its inflow. In 
this case, it is easily shown that. 


^AN * 

0 

(a) 

^TN “ 

0 

(b) 

a 

0 

(c) 

II 

r— 1 

H 

sin a 

(d) 

^TN “ 

0 for N>1 

(e) 


With this input, and first solving for the 
steady loads in order to get ((}> + Uj) , the 
results presented in figures 9 and 10 were 
obtained. It should be noted that V cos a 
must be used as the steady inflow velocity 
for the propeller at an angle of attack. 

Figure 9 presents the azimuthal variation 
of the gradient section of the thrust at the 
0.75 station for a 3-bladed propeller at a 
4.6' angle of attack. ACy, the increment in 
the thrust coefficient due to the unsteadiness 
is defined in accordance with standard term- 
inology as. 


( 18 ) 

pn D 

Three curves are sho'wn on the grrph. The ex- 
perimental curve is taken from reference 9. 

The curve labeled "predicted, 2D unsteady" is 
based on the methods just described. In B. 

figure 10, the gradient of ACj with respect 
to the dimensionless radius, x, is presented. 

The quasi-steady calculations were performed 
by simply assuming that steady aerodynamics 
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could be applied to the instantaneous section 
velocity diagraa. This is tantamount to using 
a C(k) value of^^unity in equation 12 and 
neglecting the h term in that equation. 

The simpler, quasi-steady approach results 
in predictions which are close to experimental 
results in amplitude but are shifted somewhat 
in phase. The unsteady 2D aerodynamic cal- 
culations generally have the correct phase but 
tend to overpredict the amplitudes of the 
unsteady forces. As shown, the ampli,tude of 
unsteady thrust is appreciable for this case, 
being of the order of 25% of the steady load- 
ing. 

Now consider the case of a propeller 
operating on a typical single-engine, light 
airplane, in this case, a Piper Cherokee 180. 

The Department of Aerospace Engineering at 
The Pennsylvania State University owns and 
operates this particular airplane; hence its 
choice. To begin, the cowling of the Cherokee 
was carefully measured using surveying transits 
in order to define its contour. Aft of the 
cowling, 3-view drawings were assu.ned to be 
adequate. Using this fuselage geometry, cal- 
culations of the velocity field at the pro- 
peller plane were performed for two different 
fuselage angles of attack of 8° and 2'. The 
higher angle corresponds to an airspeed of 
approximately 35.8 m/s (80 mph) while the 
lower angle is typical of cruising conditions 
of approximately 55.0 m/s (123 mph). Cal- 
culated tangential and axial velocity dis- 
tributions at the propeller plane are present- 
ed in figure 11 for the 2® angle of attack 
case. Figures 12 and 13 present the first two 
harmonic amplitudes tor both 2“ and 8® as a 
function of propeller radius. Not unexpected- 
ly, the amplitude... are seen to increase toward 
the hub. Also, observe that the axial com- 
ponent of v..iocity IS rot as dependent on a 
as the tongential component. The slight 
asynn'.-Cries about values of 0® and 180® 
showi in figure 11 result from the fact that 
the<e components are relative to the propeller 
plane which is yawed slightly with respect to 
the fuselage plane of symmetry. 

The calculation of the axial and tan- 
gential velocities in the plane of the pro- 
peller depend primarily on the shape of the B 

cowling and angle of attack. The latter re- 
sults directly in a tangential component and 
also deter^iines tlie upwash generated by the ( 

wing. Tne insensitivity of the flow field to 
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the fuselage geomt-^ry aft of the rowling is 
shown in figure lA. In one case the complete 
fuselage is modeled. In the other, the fuse- 
lage is cut off behind the cowling and re- 
placed by a smaller afterbody which can be 
defined using fewer panels. This simplification 
saves computation time and is seen to have al- 
most no effect on the velocity field except 
near the top where = 0 and the elimination of 
the windshield has tne greatest influence. 

Even here, however, the cnange in the velocity 
is only approximately 2J1. 

Using drawings of the propeller for the 
Cherokee 180 propeller supplied by the manu- 
facturer and the calculated velocity field, 
predictions have been made of the propellers 
steady and unsteady performance. Figure 15 
presents the predicted steady performance of 
the Sensenich propeller used on the Cherokee 
180. Series 16 airfoil sections were assumed 
in making these predictions since their con- 
tours can be closely matched to tne special 
sections developed and used b> Sensenich. 

These results are typical of those to be 
found in the NACA reports such as reTersnce 
10. The propeller appears to be well matched 
to the airplane in the cruise condition where 
It operates at an advance ratio of 0.75. 

Figure 16 illustrates the imsteady results 
which were obtainad for the Cherokee 180. 

Here, the first harmonic of the unsteady thrust 
gradient for angles of attack of 2® and 8“ is 
presented as a function of the radial position 
along the blade. Kote that the predicted 
amplitude of the unsteady thrust is nearly 
constant along the blade. The steady loads, 
of course, increase toward the tip but this 
is not so, at least in this example, for the 
unsteady loads. The unsteady velocity com- 
ponents apparently decrease with increasing 
radius so as to compensate for the increasing 
resultant velocities due to tor. 

PROPELLER BLADE DYNAMICS 


The unsteady aerodynamic forces just 
discussed will excite a dynamic response from 
the propeller blades. Added to this excitation 
will be, of course, any unsteadv torque im- 
pulse from the engine crankshaft. In order B. 

to predict fatigue stresses on the propeller 
blades, one must in the final analysis treat 
the complete engine-propeller d\namlc system. 

As an important part of the complete analysis. 
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a study has beer, undertaken of the structural 
dynamic behavior of the propeller as an isolated 
body. The object of this study, which has 
been accomplished, was to develop a method for 
predicting the normal modes of vibration for a 
propeller for the case where the hub is rigid- 
ly clamped. These modes are determined experi- 
mentally by propeller manufacturers to assure 
that none of the lower modes coincide with the 
exciting frequencies from the powerplant. 

However, tne testing is normally performed, 
with the propeller supported a-; a free body. 

It was hopes to model coupled in plane 
bending-out of plane bending torsion 
motion. VJhile this goal has proven to be 
elusive, a computer program has been developed 
which will predict the normal modes for coupled 
bending-bending or bending-torsion motion. 

Since propeller blades are very stiff in plane 
and torsionally the lover modes of the coupled 
bending-bending and coupled bending-torsion 
modes have approximately the same frequencies 
which are determined principally by the 
relatively soft flapwise (out of plane) bend- 
ing stiffness. Thus failure to model the 
coupled bending-bending-torsion motion is not 
considered to be too serious. 

A major accomplishment of this study has 
been the development of a program to calculate 
the shear center location for a solid cross- 
section having an airfoil-like shape. This 
program, based on a membrane analogy, and 
soap film rheory, will accept as input discrete 
points on an airfiil contour or an analytical 
definition of the coi-cour. Some results of 
this program arp ,jresented in figure 17. 

Uithin a given airfoil family, the position 
of the shear center is seen to be nearly 
constant. However, the chordwise position 
differs significantly between the NACA 4-digit 
and series 16 airfoil families. In addirion 
to the location of the shear center, this 
program also calculates the location of the 
section center of gravity, principal axes and 
ticments of inertia. All of these section 
properties are required in order to calculate 
the normal mode shapes and frequencies. 

These section properties for the Scnsenich 
propeller model 76E>'8-0-61 are presented in 

figure 18. B. VJ. McCormick 

The relationships on which the computer 
program for determining the normal nodes are 
based are too lengthy to be presented here. 

They can be found in reference 3 and depend 
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on the concept of 3 trancnission natrix. The 
original source for cucn of this analysis can 
oc found in reference 11. 

The predi .ions of ratural frequencies 
have been confimed by conparison to laeasure- 
nents provided by the oanufacturcr and by our 
own testing. The experinentai setup which was 
used is shown schenatically in figure 19. Tne 
propeller was excited as a free body by placing 
it on a soft bed nade fron two inncrtubes. In 
order to cantilever a blade, the hub was 
clanped in a Universal Testing nachine. Even 
though the nass of tnt nachine is large and 
the clacping pressure high it was neccssarv 
to weigh’ the blade opposite tne one being 
excited .r. order to avoid spurious codes. 

Soce of the experunental and calculated 
results for the Sensenich propeller are shown 
In Table 1. Generally, the frequencies which 
were obtained experimentally for the first four 
syrmetric nodes of the propeller supported 
freely confim those found by the nanufacturer . 
With the hub clanped, the predicted frequencies 
for the first nodes for coupled bending-bending 
and coupled bending-torsion are nearly equal 
and close to the frequency ncasured for the 
cantilevered blade. 

Tne predicted node shapes for the first 
three nodes for the clanped propeller are 
presented in figure 20. The prcJicte.' fre- 
ouencies foi these nodes are given in figure 
21 as a function of propeller rotational 
speed. Also included in this figure are the 
results supplied b> Stnsenich for the freely- 
supported propeller. Centrifugal stiffening 
IS seen to increase the natural frt-auencies 
of each rode slightly as the propeller speed 
increases. Superimposed on this figure are 
lines represent I'^g integer nultiples of the 
torque impulse fror the engine. Since half 
of tne cylinders of a four-cycle, horizontally- 
opposed engine fire during each revolution, 
the fundamental impulse frequency for the 
four cylinder Cnerokee ISO engine in cycles 
per second, is 

_ . no. of cvlinders , , rpm . 
p . I _ j ( --- ) Hz 


F = 


JP= 1 
60 ^ 


B. Vf. McCormick 


Tnis l..pe of diagram, ’-here ociting fre- 
quencies are super i”pcsed over natural 
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frequencies is frequently applied to rotating 
systens othei chan propellers and is knovm a^ 
a Southwell diagran. In order to avoid exces- 
sive vibration, the operational rotational 
speed of the systen, in this case the propeller 
rpis, should be selected so that none of the 
natural frequencies and exciting frequencies 
are coincident. In the case of the Cherokee 
180, the rcconaendei cruising rpa of 2500 is 
above both the first ano second nodes of either 
the cantilevered or freeiy-supported blade. 

SU>!MAR1 

A study of the aerodynanics and dvnanics 
of propellers for general aviation aircraft 
has resulted in the development of three con- 
puter progreas. The analytical bases for 
these programs can be found in references 1, 

2, and 2. In addition, user's manuals are 
currently being prepared together with pro- 
gram listings in FORTRAK. briefly these pro- 
grams 

0 calculate the velocity field in the 
propeller plane for a wing-fuselage- 
nacelle conbinatlon, 

S predict the steady and unsteady pro- 
peller aerodynamic loads, 

0 predict the natural mode shapes and 
freoLen''ics for the propeller up to 
the first four nodes. 

In the future, it is hoped to combine the 
programs in order to predict the forces 
exciting airframe vibration. A normal nedes 
approach will probably be tried to accomplish 
this challenging task. In addition, it is 
also planned to obtain additional experinenfal 
confirmation of the existing programs. 
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Table 1 


NATURAL FRCOUENCIES OF A NON-ROTATING SEIISENICH PROPELLER 



Free Beam S 

BffiRBRKIIBlfBSiil 

Cantilever Beam | 




Predicted Values I 

Experimental Data 

Mode 

Number 

Scnnenlch 

Data 

Experimental 

Data 

Coupled 

DendlnR-Bcndin 


Coupled 

BendlnR-Toralonal 


1 

56.8 

58.2 

51.7 


51.4 

50.0 

2 

170.0 

175.0 

161.5 


168.1 

139.0 

3 

3A2.9 

3A7.0 

312.7 


372.3 

316.0 

A 

A86.6 

1 

431.0 

384.8 


664.0 

434.0 










FIG. 1 PANELING OF BODY SURFACE 












PROPELLER NACA 5868-9 
CLARK Y SECTION 
NO OF BLADES = 2.3 & 4. 
BLADE ANGLE AT 3M STATION 
ACTIVITY FACTOR « 78.9 



L.Z 1.6 1.2 1.6 1.4 

ADVANCE RATIO J 

FIG. 4 COMPARISON OF PREDICTED AM) f'EASURED PROPELLER 
EFFICIENCY FOR DIFFERENT NUMBER OF BLADES 





POWER COEFFICIENT C 


II 


tCL 

PROPELLER NACA 10- (3H05) -045 
NACA 16-SERIES SECTION 

NO OF BLADES =2 ... 

BLADE ANGLE AT 3/4 STATION =20, 25 & 35 
ACTIVITY FACTOR = 131.5 
ROTATIONAL SPEED = 2000 RPM 



FIG. 5 COMPARISON OF PREDICTED AM) MEASURED 
POi'/ER COEFFICIENT AT 2000 R.P.M, 


I 

J 


MACH NUMBER M S M, 





FIG. 6 VIEV/ OF PROPELLER LOOKING 
IN DIRECTION OF FLIGHT 



th 


FIG. 7 H 


HARi'-IONIC OF INFIJO.-/ VELOCITY FOR BLADE SECTION 



Q 




CHANGE IN SECTION THRUST 
COEFFICIENT 6RADI ENT ACT/dX 


propeller NACA 10- (3) (08) -03 
NACA 16-SERlES SECTION 
NO OF BLADES = 3 
BLADE ANGLE AT 3/ A STATION = 30* 
ACTIVITY FACTOR = 86.2 
PROPELLER ANGLE OF ATTACK = A. 55* 
ROTATIONAL SPEED = 1350 RPM 
ADVANCE RATIO = 1.2 



FIG, 9 CHANGE IN BLADE SECTION THRUST COEFFICIENT GRADIETiT WITH 
AZIfWIAL POSITION OF BLADE AT THE 3/A RADIAL STATION 


PERCENTAGE OF 
STEADY STATE VALUE 



AMPLITUDE OF UNSTEADY SECTION 
THRUST COEFFICIENT GRADIENT d( C_)dX 



PROPELLER NACA 10- (3) (08) -03 

NACA 16-SERIES SECTION 

NO OF BLADES = 3 

BLADE ANGLE AT 3/4 STATION = 30 

ACTIVITY FACTOR =86.? 

PROPELLER ANGLE OF ATTACK = 4.55 
ROTATIONAL SPEED = 1350 RPM 



FIG. 10 DISTRIBUTION OF Af^LITUDE OF UNSTEADY SECTION 
THRUST COEFFICIENT GRADIEITT 










AMPLITUDE OF HARMONICS OF AXIAL VELOCITY (V-/v 



FIG. 12 HAR’IONIC CaiTENT 0- THE AXIAL VELOCITY DISTRIBUTION 
IN THE PROPELLER PLA*.E CPEP<ATII.'G IN THE FLO,/ FIELD 
OF A OPER CHEROKEE 160 





AXIAL FREE STREAM 



COWLING ANGLE-OF-ATTACK =2.0° 



FIG, 14 CQ’PARISON OF AZIfVIFAL DISTRIBUTIONS OF AXIAL 
VELOCITY IN THE PRC?ELL£R PLA'-'E OF A PIPER 
ChEROr^E 180 PREDICTED USING 1V.0 DIFFERErfT 
BODY GEaiETRY MODELS AFT OF THE CWLING 




PROPELLER SENSENICH 76EM8-0-61 
EQUIVALENT NACA 16-SERlES SECTION 
NO OF BLADES = 2 

BLADE ANGLE AT 3/4 STATION =18.5 
ACTIVITY FACTOR = 94. A 
ROTATIONAL SPEED 2700 RPM 





FIG. 15 PREDICTED STEADY STATE PERFORfWJCE OF THE 
SENSBJICH PROPELLER 76E/'iS-0-61 


EFFICIENCY 



AMPLITUDE OF UNSTEADY SECTION 
THRUST COEFFICIENT GRADIENT d(/^Cy)/dX 


St 


PROPELLER SENSENICH 76EM8-0-61 
EQUIVALENT NACA 16-SERIES SECTION 
NO OF BLADES = 2 

BLADE ANGLE AT 3/4 STATION =18.5 
ACTIVITY FACTOR =94.4 
ROTATIONAL SPEED = 1900 RPM 



0.0 0.4 0.6 0.8 1.0 

NONDIMENSIONAL RADIAL STATIONyX 


FIG. 16 DISTRIBUTION OF AMPLITUDE OF UNSTEADY SECTION THRUST COEFFICIENT 
GRADIENT FOR THE PROPELLER IN THE FLOW FIELD OF THE CHEROKEE 180 




CHORDWISE POSITION OF THE SHEAR CENTER. AND CENTER OF GRAVITY 

EXPRESSED AS A FRACTIOn'oF CHORD 
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FIG, 17 SHEAR CEI.TER LOCATION FOR NACA 4-DIGIT SERIES 
VERSUS NACA 16 SERIES AIRFOILS 


r = CENTER OF GRAVITY LOCATION 

Xe r = SHEAR CENTER LOCATION 

O t w t 


NACA 16 OXX SERIES SYMMETRICAL AIRFOILS 


NACA 4 - DIGIT SERIES SYMMETRICAL AIRFOILS 
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MAXIMUM THICKNESS RATIO, (t/c). 
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FIG. 18 I-TATERIAL AND GEQ"iE:TRICAL PROPERTIES OF A SENSENICH PROPELLER BLADE 
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AREA MOMENT OF INERTIA OF THE BLADE SECTION 
ABOUT MAJOR PRINCIPAL CENTROIDAL AXIS> 

X 10 (n^) 

Blade angle relative to the chord line> 

6 (deg) 


distance between mass 
AND elastic axis POSSITIVE 
WHEN mass lies AHEAD OF 
THE SHEAR CENTER^ 

E X lo"^ (m) 
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AREA MOMENT OF INERTIA OF THE BLADE SECTION MASS PER UNIT LENGTH^ 

ABOUT MINOR PRINCIPAL CENTROIDAL AXIS^ M (kG^M) 

. • , --6 / 
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FIG. 19 


SCHEmilC OF TEST SETUP FOR DETERMINING 
PROPELLER BLADE NATURAL FREQUENCIES 







FLAPWISE MODE SHAPE^ W 



FRACTIONAL RADIUS^ X 

FIG. 20 PREDICTED MODE SHAPES OF COUPLED FUU’.'/ISE 

BENDING AND CHORWISE BENDING VIBRATIONS FOR 
A CLAJ-'PED SENSED I CH PROPELLER 


NATURAL FREQUENCY/ 



FIGi 21 PROPELLER NATURAL FREQUBJCIES AND ENGIfiE EXCITING FREQUENCIES 
AS A FUNCTION OF PROPELLER SPEED 
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